Quantum magnetooscillations in the ac conductivity of disordered graphene 
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The dynamic conductivity a(oj) of graphene in the presence of diagonal white noise disorder and 
quantizing magnetic field B is calculated. We obtain analytic expressions for o-{ui) in various para- 
metric regimes ranging from the quasiclassical Drude limit corresponding to strongly overlapping 
Landau levels (LLs) to the extreme quantum limit where the conductivity is determined by the 
optical selection rules of the clean graphene. The nonequidistant LL spectrum of graphene renders 
its transport characteristics quantitatively different from conventional 2D electron systems with 
parabolic spectrum. Since the magnetooscillations in the semiclassical density of states are anhar- 
monic and are described by a quasi-continuum of cyclotron frequencies, both the ac Shubnikov-de 
Haas oscillations and the quantum corrections to o-(u) that survive to higher temperatures manifest 
a slow beating on top of fast oscillations with the local energy-dependent cyclotron frequency. Both 
types of quantum oscillations possess nodes whose index scales as uj 2 . In the quantum regime of 
separated LLs, we study both the cyclotron resonance transitions, which have a rich spectrum due 
to the nonequidistant spectrum of LLs, and disorder-induced transitions which violate the clean 
selection rules of graphene. We identify the strongest disorder-induced transitions in recent magne- 
totransmission experiments. We also compare the temperature- and chemical potential-dependence 
of o(uj) in various frequency ranges from the dc limit allowing intra-LL transition only to the uni- 
versal high-frequency limit where the Landau quantization provides a small B-dependent correction 
to the universal value of the interband conductivity a — e 2 /Ah of the clean graphene. 

PACS numbers: 72.80.Vp, 73.43.Qt, 78.67.Wj 



I. INTRODUCTION 

Since its discovery in 2004pJ the two dimensional carbon 
allotrope - graphene - is attracting outstanding interest 
in the condensed matter community. It has been verified 
that carriers in graphene show a linear dispersion relation 
with Fermi velocity vq « 10 6 m/s and are governed by the 
massless Dirac equation!^ Due to the carriers' Dirac na- 
ture, graphene shows remarkable properties. The clean 
density of states (DOS) is linear in energy and vanishes at 
the Dirac point, vv{e) = H/27™ 2 (h = 1). The absence of 
scales at the Dirac point gives rise to a universal dc con- 
ductivity of the order of the conductance quantum. The 
universal high-frequency conductivity a = e 2 /Ah yields 
the constant absorption coefficient of 2.3% which makes 
graphene attractive for broadband optical applications. 4 
The nontrivial topology of graphene leads to the char- 
acteristic half-integer quantum Hall (QH) effect and to 
the nonequidistant Landau level (LL) spectrum includ- 
ing the unusual zeroth LL at e = 0. 6 ' Due to the large 
cyclotron energy lu c = v W2eBjc compared to conven- 
tional 2D electronic systems (2DES) the QH effect can 
be observed up to room temperature.^ Recent rapid de- 
velopments demonstrate great potential of graphene in 
optoelectronics PEl From this perspective, the applica- 
tion of a quantizing magnetic field creates a suitable en- 
vironment for applications in e.g. laser physics!^ The 
quantization introduces a tunable energy scale in the oth- 
erwise scale-free graphene, while the nonequidistant LL 



spectrum makes it highly selective in the frequency do- 
main. 

While experimental data on quantum magnetooscil- 
lations in graphene is limited (although growing) DUEI 
there are comprehensive studies of related phenomena for 
semiconductor 2DES where the effects both near and far 
from equilibrium have been extensively studied! 18 ^ 19 ! In 
the linear response regime, quantum magnetooscilla tions 
in the ac conductivity were theoretically predicte d 20 * 21 ^ 
and observed^ already in the seventies. In Ref. [53] 
the theory was generalized to high mobility 2DES with 
smooth disorder potential; the findings are also confirmed 
in recent experiment! 2 ^ 

This work is (i) motivated by experimental and tech- 
nical advances in the graphene research and (ii) gen- 
eralizes the theory on quantum magnetooscillations in 
the conventional 2DES mentioned above. We study 
quantum oscillations in Landau quantized graphene in 
the presence of disorder. Previous theoretical works 
on the magnetoconductivity in graphene accounted for 
the disorder in terms of a phenomenological broaden- 
ing in the single-particle s pectr urrl 2 ^ 2 ^ or focused on 
the dc magnetoconductivityPSM! The optical conductiv- 
ity has also been studied for vacancies within the T- 
matrix approximation.^ Here we perform a systematic 
calculation of the dynamic magnetoconductivity within 
the self-consistent Born approximation (SCBA). 

The paper is organized as follows. In Sec. [H] we out- 
line the model of the SCBA in graphene and present re- 
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suits on the spectrum of disordered graphene in a mag- 
netic field, in the semiclassical as well as in the quantum 
regime. Further details on the calculation of the density 
of states can be found in Appendix[A| Scction |ITI| presents 
the formalism for the calculation of the dynamic conduc- 
tivity. The following Sees. |IV| and [V] are devoted to the 
dynamic conductivity in the regimes of strongly over- 
lapping and well-separated LLs respectively. In Sec. |VI| 
we calculate the high-frequency conductivity of a mod- 
erately disordered graphene. The summary of results is 
presented in Sec. |VII[ Details on the calculation regard- 
ing Sees. |IV| and [V| can be found in Appendix [B] 



II. LANDAU LEVEL SPECTRUM WITHIN 
SCBA 

In the following we assume that disorder does not mix 
the two valleys in graphene and therefore calculate all 
quantities per spin and valley. Correspondingly, to get 
the full conductivity of graphene one should multiply the 
results for conductivity below by the degeneracy factor 
of 4. In the presence of a constant magnetic held in z- 
direction the electrons in a single valley of clean graphene 
are described by the 2D Dirac equation 



H 



v a ■ 



p— eA(r) , A(f) = —yBx 



(1) 



Here we have chosen the Landau gauge for the vector 
potential A. The vectors r — (x,y) T , x — (1,0) T , and 
<? = {pxi^y) denote the Pauli matrices. The positions of 
the Landau levels (LLs) in clean graphene are given by 



E n = sign(n)w< 



n G Z 



(2) 



The optical selection rules of the clean graphene allow 
transitions from the LL n to m if 



±1 



(3) 



therefore enabling both intra- [sign(n) = sign(m)] and 
interband [sign(n) =^ sign(m)] transitions. Since the LLs 
move closer at higher energy, it is convenient to introduce 
a local cyclotron frequency 



E\ 



\n\ + l 



E\ 



!»! 



(4) 



that is, the distance between neighboring LLs. In high 
LLs, n 1 (s ^> w c ), it approaches the quasiclassical 
cyclotron frequency of a massless particle 



Within the SCBAj^EU the self-energy is given by 



iq-r 



(7) 



Weak disorder is characterized by the disorder correlator 
W{q) which in the case of graphene is generally a fourth 
rank tensor in valley and sublattice space. Here we con- 
sider short-range impurities which do scatter between the 
sublattices but do not produce any intervalley scattering. 
In a given valley, this gives diagonal white noise disorder 
with the correlator 



(8) 



characterized by a single parameter a -C 1. For the di- 
agonal disorder, the self-energy is independent of the LL 
index and becomes diagonal in the sublattice space; still, 
it carries an asymmetry between the sublattices a and 6, 



£ = diag(S a ,E b 



(9) 



The asymmetry is due to the fact that - in a given valley 
of the clean graphene - the wave function of the zeroth LL 
resides in one sublattice only. In what follows, we choose 
the valley such that the wave function of the clean zeroth 
LL resides in the sublattice b. The SCBA equation ^ 
for disorder with the correlator (|8| acquires the form 



Ea(6) = 



E 

1=1(0) 



-6(a) 



(e-E Q )(e-S fc )- Wc 2 7 



(10) 



Away from the zeroth LL the difference between the two 
self-energy components is negligible, E Q ~ £;,, [see dis- 
cussion under Eq. ( 16 ) and Fig. [I] yielding 



2 +°° 



(e-E) 



(11) 



Apart from specifics of the zeroth LL due to its pro- 
nounced sublattice asymmetry, one can consider two lim- 
iting cases we address separately below: (i) clean, or 
quantum limit, when the disorder-broadened LLs remain 
well separated and (ii) dirty, or classical limit when LLs 
strongly overlap almost restoring the linear slope of the 
DOS at B = 0. Unlike conventional 2D systems with 
parabolic spectrum, in graphene the two cases (i) and 
(ii) frequently coexist: LLs, well separated near the Dirac 
point, start to overlap at higher energies where the local 
cyclotron frequency (|5| strongly reduces compared to lo c . 



,loc 



2kl 



(5) 



A. Separated Landau levels 



The disorder is included into the self-energy E which 
enters the impurity averaged electronic Green's function 



G=(e-H 



(G) 



We start with the limit of well separated LLs. In this 
case, the main contribution to the self-energy at energy 
e comes from the states in the nearest LL of the clean 
graphene to which we assign the integer number N closest 
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to e 2 /uj 2 . LLs with index n ^ N contribute to logarith- 
mic energy renormalization as detailed below. 

For N 0, the sublattice asymmetry can be neglected, 
and the solution to Eq. (11) for the retarded self-energy 
reads [the advanced self-energy T, A = (£"**)*] 



e-E H = 



E 



N 



where the width of the iVth LL 



fauj c /Z{E N ) , N ^ 



N^O, (12) 



(13) 



and the condition of applicability is Tn <C i^n- Apart 
from the usual renormalization of energy by the factor of 
2 inside the LLsj 21 * 31 ^ in graphene the energy gets addi- 
tionally rcnormalized according to 



i = e/Z(E N ). 
Here the renormalization constant is 

Z(e) = l-aln(A c /|e|), 



(14) 



(15) 



and A c is the high energy cut-off of the order of the 
band width. Within the SCBA, the additional logarith- 
mic correction describes the influence of states in dis- 
tant LLs [n ^ N in Eq. ([10|], similar to renormaliza- 
tion group (RG) corrections P For the type of disorder 
we are investigating here, the SCBA corrections with- 
out magnetic field are known to be in quantitative ac- 
cordance with RG calculations.^ Apart from the addi- 
tional logarithmic renormalization and from the non- 
equidistant spectrum E^ of clean LLs, Eq. ( 12 ) repro- 



duces the well-known semicircular law o btained by Ando 
for 2DES with parabolic spectrum)2il21lFrom the require- 
ment aln(A c /|e|) <C 1, we obtain that the results are jus- 
tified for energies above the exponentially small energy 
scale Ace' 1 / 01 . 

In the vicinity of the zeroth LL, N = 0, one needs 
to take into account the explicit sublattice asymmetry. 
Equation ( 10 ) yields the self-energies^ 



e i 



2 2 



(16) 



The width of the zeroth LL in both sublattices is 



(18) 



the factor V2 comes from different degeneracy of the ze- 
roth LL - the total number of states A/tot = L 2 /2irl 2 B is 
the same in all LLs, but in the zeroth LL all these states 
reside in one sublattice. We observe that only a small 
amount of the spectral weight oc 1 — Z(E\) oc a <C 1 
is scattered into the afore empty sublattice a. At the 
same time, the strong renormalization of energy by the 
factor of 2 related to the high degeneracy of clean LLs in 
Eqs. (fl2]) and (Il6] is absent in Eq. ([17}. Only the log- 



arithmic renormalization due to distant LLs with n 7^ 
remains. 



B. Overlapping Landau levels 

In view of the condition a < 1 of weak disorder, the 
regime of strongly overlapping LLs is realized at high 
energies e > e ov ~ ui c /y/a ^> lu c . Indeed, the LL number 
-ZV ov where LLs start to overlap is given by the condition 
Tn ~ oj n which gives [see Eqs. ©, @, and ffify] 



a 1 > 1. 



(19) 



Accordingly, we assume T, a j, ~ S and use the Poisson 



summation formula to rewrite Eq. (Ill into a rapidly 



convergent sum in the Fourier space [see Eq. (Al)]. The 



small parameter that controls such an expansion is the 
Dingle or coherence factor 



A 



(20) 



which is the analog of the Dingle factor 5 = c^ n ^ cTq in 
2DES with parabolic spectrum and describes the smear- 
ing of quantum oscillations in the disordered system. To 
zeroth order in A, we obtain a broadening of the single- 
particle states at B = 0, 



lA-s-0 



= -i/2r,(e). 



(21) 



with the energy-dependent quantum scattering time 
given by 



T q {e) 



Z(e) 



COT £ 



(22) 



The renormalization constant Z(e) which defines the 
renormalized energy e — ejZ(e) is given by Eq. (15 1. 



To first order in A 1, the self-energy acquires the form 



2r g (e) 



2Ae i27rel£l/ ^ 2 + 1 



(23) 



C. Density of states 

The DOS per spin and valley v(e) = — Im tr G R (f, r, e)/ir 
in the case of white noise disorder is given by 



V {e) 



1 



2n 2 v 2 a 



Im trS R 



(24) 



The DOS in the sublattice b obtained numerically from 
10 1 is plotted in Fig. [T] as the thin 



the SCBA equation ( 
line. Separated LLs s 



low a semicircle DOS, see Eqs. ( 12 ) 



and (16). The dashed line for the DOS in the sublattice 
a shows that (i) in the zeroth LL the DOS transferred 
into the sublattice a is small and vanishes if disorder is 



turned off, see Eq. ( 17) and (ii) at higher energy the DOS 
(self-energy) in both sublattices is approximately equal. 
Finally, the thick line, calculated according to Eqs. ( [23] ) 
and ( 24 ) , illustrates the limit of strongly overlapping LLs 
corresponding to large e >• e v (which gives A <?C 1). 
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FIG. 1. The DOS in disordered graphene. The thin line is the 
DOS in t he su blattice b obtained numerically from the SCBA 
equation (10 1 for a = 0.01. For the first two levels we show 
the DOS in the sublattice a as the dashed line. Already in 
the second LL there is no noticeable difference between the 
a and b sublattice. The thick curve shows our result for the 
DOS according to Eqs. (|23[) and (241. It is in good agreement 



with the numerical results for high energies. 



D. Semiclassical regime 

Here we address the semiclassical regime of high LLs 
in graphene, N ^> 1, additionally assuming that these 
levels strongly overlap, e e ov ~ u> c /y/a. The aim 
is to provide a link to the the results in the corre- 
spondent parabolic band with effective energy-dependent 
mass such that the local characteristics of graphene and 
conventional 2DES with such parabolic spectrum at high 
LLs are identical. 

In the limit A 0, Eqs. ([211 and pi} yield the B = 



renormalization of the clean DOS vq by disorder 



2irv$ 



v = Z(e)u(e)\x^o 



According to Eq. ( 23 ) , the correction oc A reads 



v{e) — i>o — 2i>qX cos 



2ne\e\ 



(25) 



(26) 



The energy renormalization induced by the Landau quan- 
tization shows similar (but phase-shifted by ir/2) oscilla- 
tions 



AReE(e) = 



A 



2ne\e\ 



(27) 



Apart from (i) the large-scale logarithmic renormaliza- 
tion ( 25 ) which is specific for the linear spectrum of 



graphene and (ii) a phase shift related to a nontriv- 
ial Berry's phase of graphene the local characteristics 
of graphene and 2DES with parabolic spectrum at high 
LLs are identical as expected. Namely, quasiparticles in 
graphene at energy e close to e* 3> lo c behave equivalent 



to massive electrons with mass, energy, and cyclotron 
energy given by 



6^ = 8/2-^/2, 

c4 eff ) = eB/m eS c, 



(28) 
(29) 
(30) 



such that the local velocity in the parabolic band 

W e il = V^off/™eff coincides with vq. The term —uJf^/2 
accounts for the shift of LLs due to vacuum fluctuations 
in the parabolic band which is absent due to the nontriv- 
ial Berry's phase in graphene, see also Eq. (33) below. 

The effective cyclotron frequency cuc"^ coincides with 
the renormalized local cyclotron frequency in graphene, 
Q l ° c (e*) = uji/2\e*\ = Z(e*)oj l ° c (e*). If one addition- 
ally introduces the renormalized f q (e) = T q (e)/Z(e) = 
l/Trajfl, equation ( |26| transforms into conventional ex- 
pression for parabolic band, specifically, 



e* 


TOeff 


2ttVq 




A = c - 2m ^ 




2Tre\e\ 


27re e ff 




,,(eff) + 



-7r/tD( ef£ )f a (e) 



(31) 
(32) 
(33) 



III. DYNAMIC CONDUCTIVITY 

Below we calculate the dynamic conductivity of disor- 
dered graphene in the presence of a quantizing magnetic 
field. We use the Kubo formula for the real part of the 
diagonal conductivity 



d£ h - /, 



■in 



K(e,e + uj) 



(34) 



to calculate the magnetoconductivity in the linear re- 
sponse. Here f £ is the equilibrium Fermi-Dirac distribu- 
tion function, and the conductivity kernel is given by 



K{s 1 ,e 2 ) = -tr 



From the Hamiltonian ([l]), the current density operator 



J 



e dH 
L dk 



evo 
L 



a, 



(36) 



where L is the length of the system. Due to the linear 
spectrum the current operator does not depend on the 
magnetic field. Since the two valleys are decoupled, we 
calculate the conductivity per spin and valley. 

The effect of the disorder averaging in Eq. (35) is 



twofold: (i) the bare Green's function is replaced by the 
impurity averaged Green's function (J6|, (ii) the summa- 
tion of the diagrams 11 'c) leads to vertex corrections 
to the current operator (36). In 2DES with parabolic 



spectrum and white noise disorder, the vertex corrections 



5 



are absent. By contrast, in graphcne the vertex correc- 
tions are present for the diagonal white noise disorder as 
well. They originate from the nontrivial Berry's phase of 
Dirac fermions. We present details on the calculation of 
the conductivity including vertex corrections in App. [B] 
As expected, we find that in the quasiclassical regime 
of strongly overlapping LLs the vertex corrections give 
rise to the replacement of r q by the transport scattering 
time Tt r = 2r q in the Drude part of the conductivity, see 



Eq. (41) below, while the quantum time appears only in 



the quantum corrections related to the Landau quanti- 
zation. In particular, r q enters the Dingle factor A, see 



Eq. (20) 



It follows from Eq. ( 23 ) that to zeroth order in A only 



the RA-sector of Eq. ( 40 ) contributes to the conductivity 
for the intraband processes. For uj = \e± — e 2 \ <S £f the 
kernel ( 39 ) has a Drude form with the energy-dependent 
broadening ( 22 ) and local cyclotron frequency ^ . If the 
temperature T -C ep, these parameters can be evaluated 
on-shell. The conductivity <td = <Jd,+ + a D,- m this 
regime reads 



<7n + (UJ) = -, 

' 4tt ( w ± uj}° c ) 2 + t-.- 



(41) 



IV. OVERLAPPING LANDAU LEVELS 



Due to the vertex corrections the transport scattering 
time r tr = 2r q replaces the quantum scattering time 
as explained in App. [B] The Drude weight V for Dirac 
fermions is 



In this subsection we consider the case of highly doped 
graphene with the Fermi energy ep 3> lj c ,lj. Therefore, 
only intraband processes are possible and one can neglect 
the sublattice asymmetry in the self-energy. It follows 
that the conductivity in this regime should not change 
in the presence of intervalley scattering. Using the conti- 
nuity equation, the conductivity kernel ( 35 1 is expressed 



in terms of density-density correlators II. Their general 
form is given in Eq. (B14), which simplifies to 



ttRR(RA) 



n=0 



Gn+l-( £ l) + ^n+l,+ ( £ l) 



G R n ^(e 2 ) + G^\e 2 ) 



(37) 



if the sublattice asymmetry is absent. Here we intro- 
duced the chiral Green's functions 



G 



R{A),~ 



n,± 



1 



e - lmZ R ( A ) =F \E„ 



(38) 



In the following we write e for e for brevity. In terms of 
the Il-correlators, the kernel (35) acquires the form 



if(£ 1 ,£ 2 ) = e 2 Re[([nf i ^ 2 ]- 1 

_! ' (39) 

-([n£? ea ] X -<wr) +{£i^£ 2 } 

As previously, we use the Poisson formula to rewrite the 
sums occurring in Eq. 



the Fourier space. Wit 



37 1 as rapidly convergent sums in 



1 the self-energies for overlapping 



LLs from Eq. ( 23 1 , we obtain 

2(ei-££)(e 2 -E 



62 J 



( Bl _Eg)"_( Ba _E 



|sign( ei ) r^E* - sign(£ 2 )r g , e2 Ef 2 ^)| 



(40) 



We organize the following analysis in orders of A. Details 
of the calculation are presented in Appendix |B| 



V = 



e 2 \e F \ 



(42) 



which produces the standard Drude weight T> = e 2 e c s 
in the correspondent parabolic band, see Eq. (29). The 



Drude part of the conductivity in graphene was also ob- 
tained in Refs. [3U and 133] Without magnetic field the 
Drude form has been confirmed experimentallj^ though 
deviations (possibly due to interactions) were also ob- 
served. 

On top of the semiclassical Drude conductivity quan- 
tum oscillations are superimposed. There are two major 
and competing damping mechanisms present. Finite T 
leads to the thermal damping of the Shubnikov-de Haas 
(SdH) oscillations in the ac- and dc-response. Scattering 
off disorder also smears quantum oscillations and this is 
captured by the coherence factor A. The latter mech- 
anism is dominant for T below the Dingle temperature 
T D = 1/2ttt ? P 

The leading order quantum corrections which are 
strongly damped by finite temperature describe SdH os- 
cillations in the dynamic and dc conductivity. The cor- 
responding contribution to the kernel K = K + + 
Eq. (35), reads 



1 

a± 



2a| 

a3T 



Vq SC {£2) ^osc(£l) 

^0(^2) v {si) 

ARe£(£ 2 ) AReS(£i) 



(43) 



1/t«j(£ 2 ) l/r q (ei) 



where we used the abbreviation a± — T tI (uj ± uj 1 ° c ). 
The first term in Eq. ( 43 1 is produced by the quantum 



correction to the DOS, Eq. (26), while the second term 



accounts for the energy renormalization, Eq. (27). The 



first term is dominant away from the cyclotron resonance, 
see Fig. [2j In this case a± 3> 1, and the kernel (43) 
acquires the form 



goscCgjQ ^ OSC (gl) 

vq{s 2 ) i>o(ei) 



(44) 
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FIG. 2. The prefactors from Eq. (|43l 
a-/(a 2 _ + cj 2 /e 2 f ), (c) 24/(4 +T), 



(a) 24/(4 + 1), (b) 
(d) a+/(4 + w 2 /4), 
where a± = r tr (aj ± 4° c )- All curves are calculated for 
ef/uj c = 6 and a — 0.01. Away from the cyclotron resonance, 
a± 3> 1, the subleading effects due to the energy renormal- 
ization (thin lines) can be ignored. 



This form of the kernel is expected from a golden 
rule consideration where the conductivity (in the clas- 
sically strong field, a± ^> 1) is determined by the 
product of initial and final DOS v{e)v(e + uS). To 
first order in the coherence factor A, this yields 
v(e)u(e + u>)-vl~ v [v osc (s + u) + v osc (e)] oc A. 

For small T,w<ej?we evaluate the smooth functions 
&d.± & n d A on-shell but keep the energy dependence in 
the rapidly oscillating parts while calculating the con- 
ductivity (34). The asymptotics of the resulting Fresnel 



integrals provides a reasonably simple expression for the 
correction to the Drude conductivity away from the cy- 
clotron resonance at T = 0, 



<j±\uj) = 4ctd,±Acos 



2tt(£ 



sm 



2lTUJ /uj 



loc 



(45) 



Equation ( 45 ) describes SdH oscillations in the dynamic 
conductivity illustrated in Fig. [3] 

The fast harmonic oscillations with the local cyclotron 
frequency Wc° c ' are similar to those known for systems 
with a parabolic spectrum except for the absence of the 
vacuum shift due to the Berry's phase of Dirac fermions 
and hence a shift in the zero of energy to ep/2 — u>2 . 

However, the magnetooscillations in Fig. [3] further 
show a slow modulation on the scale of the cyclotron 
frequency uj c at the Dirac point. This beating in the 
quantum oscillations is due to the difference of the cy- 
clotron frequencies at the initial and final state of the 
optical transition. Remarkably, the quantum oscillations 
show nodes stemming from destructive interference of the 
density of states in the initial and final state. If the Fermi 
energy ep is situated in the center of a LL, the node oc- 
curs at 



UJ 



v / 2fcTl 



k £ N 



(46) 



FIG. 3. The SdH oscillations in the dynamic conductivity 
(solid line) and the envelope function showing slow modu- 
lation due to nonequidistant spectrum of LLs (dashed line) 
calculated according to Eq. (451 for ef/uJc ~ 6. 



Detuning of the Fermi energy from the center 
of the LL shifts the frequency nodes such that 
e% + uj 2 = uJ 2 (n/2 + 1/4), n £ N holds. Note that al- 



though we assumed U) <?C £p the frequencies in Eq. ( 46 ) 



are well within the range of our approximation for not 
too large k in view of Ep 3> ui c - For the values of uj 
in Eq. ( 46 ) the DOS oscillations acquire a relative phase 
shift of approximately (2k + l)w between e and e + uj 
due to the e-variation of the period which leads to the 
destructive interference. 

At T > Tps-, the temperature smearing dominates over 
the quantum mechanical smearing. However, there are 
additional quantum oscillations oc A 2 that survive higher 
temperatures and become exponentially larger than SdH 
oscillations at T ^> Tp> ■ Th ey ar e well known for systems 
with a parabolic spectru m 1 21 ! 23 ! and originate from the 
terms in the i?^4-sector in Eq. ( 40 ) which do not oscillate 



with mean energy e\ +£2 but do oscillate with the energy 
difference u) = |ei— £2!- Such w-oscillations are insensitive 
to the position of the chemical potential with respect to 
LLs. Therefore, averaging over the temperature window 
according to Eq. ( [34] ) does not lead to additional damping 
in contrast to SdH oscillations. The relevant contribution 
at T > T D is 



ai 9 V) = 2a D ,± A 2 



2a ± 



1 



2ttuj 



al + l 



sm 



UJ, 

2ttu 



loc 



2ituj 



uj 1 c oc J 2ttuj 2 juj 2 c 



(47) 



Away from the cyclotron resonance, aj. » 1, the quan- 
tum correction (47) acquires a simpler form, 



2 <7n ± A 2 



2nuj 2 /u> 2 



cos ■ 



2ttuj 

,loc 



(48) 



The conductivity at high temperature T 3> Tp> is plot- 
ted in Fig. [4] Apart from the integer cyclotron resonance 
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FIG. 4. Magnetoconductivity (Td,±[u) + er±(u)) of graphene 



at high temperatures [normalized to the Drude value ( 41 1] 
calculated according to (48 1 for Ef/u c = 6 



Vertical lines 

mark the position of integer harmonics of the cyclotron reso- 
nance. 



harmonics uj — nw c oc , n € N, that have an analog in sys- 
tems with a parabolic spectrum, we encounter again an 
additional modulation of the quantum oscillations due to 
the presence of multiple cyclotron frequencies. Since the 
temperature-stable quantum corrections oc A 2 are insen- 
sitive to the position of the Fermi energy wi th re spect to 
LLs, the positions of the nodes, uj/u c = y/n/2, n e N, 
are therefore solely determined by the probing frequency 



V. SEPARATED LANDAU LEVELS 

We now address the quantum regime of well resolved LLs, 
which involves LLs with numbers \n\ < N ov , see Eq. ( 19 1. 
In the following, Lk denotes the K-th LL. 

Our results for the dynamic conductivity at T = are 



illustrated in Fig. [5] As long as LLs are separated, the 
conductivity is the sum of contributions from individual 
transitions between Lk and Lm- Despite LLs are sepa- 
rated, the corresponding peaks in a(uj) overlap. Indeed, 
some transitions are degenerate (for instance Lq — > L\ 
and L_i — > L in Fig. [5]). In other cases, the excita- 
tion energies for different transitions may become close 
to each other due to the nonequidistant spectrum of LLs. 



From the Kubo formula (34) we find that to leading 



order in T/uj c the conductivity kernel is proportional to 
the density of the initial and final states. In this section, 
we neglect the renormalization of the energy described by 
Z(e), Eq. (15), assuming that the contribution of states 
with energies |e| < Ace -1 / 
Z{e) = 1 in Eqs. Q-Q, which gives ImSf 



is negligible. We thus put 
in 



Eq. \17\. Using Eq. (24), we obtain the partial DOS 



l 2 B n 2 T 2 



Rey/T 2 - (e - E K ) 2 , r = y/ai 



(49) 



in either sublattice a or b for K ^ 0, and 



v K =o(e) 



1 



z|7r 2 r 2 



Rev^r 2 - {ef 



(50) 



residing in sublattice b for K = 0. The total DOS, in- 
cluding both valleys, spin components and sublattices is 
given by 8vk^o f° r K ^ and by Avk=q for K = 0. The 
width of the individual resonances is then 2(rjy + ^k), 
determined by the width of the DOS in Lk and Lm- 
In the following we distinguish the cyclotron resonance 
peaks, \\K\ — \M\ \ — 1, from the disorder-induced peaks, 
\\K\ — \M\\ 1. The latter vanishes if the disorder is 
switched off, while the former survive as they respect the 
clean selection rules in graphene. 

Using u c S> T for the II-correlators (40), we cast the 



conductivity ( 34 ) in the form 



a(u)=a £ P^Mj^J^/i.T), (51) 



K.M 



where <j = e 2 /47r 2 and ck — 1 + Sk,o, see App. |B|and[c] 
for details. The coefficients 



V(K,M) = 1, 



\K\-\M\\ = 1, 



(52) 



for the cyclotron resonance transitions; for the disorder- 
induced transitions not involving the zeroth LL (\\K\ — 
\M\\ ^ 1, K ^ and M £ 0), 



V(K, M) 



T*(E M + E K ) 



2(E 2 



K 



E 2 

M 



+ {M++K}; (53) 



finally, for the disorder-induced transitions involving the 
zeroth LL (| \K\ - \M\ \ ^ 1, K = or M = 0) 



V{K,M) = 



r 2 (E 2 K + u 2 ) 

2{E 2 M -E 2 K +^) 2 



{M -H> K} . (54) 



The function J- km in Eq. (51) describes the shape of 



the peaks in the conductivity; its general form is given 
in App. [C] see Eq. (C5). In the case of the disorder- 
induced transitions ■ \M || 7^ 1, the vertex corrections 
are negligible: the conductivity kernel can be calculated 
using the bare polarization bubble and depends on energy 
only via the product vk{£)vm[£ + Correspondingly, 
for \\K\ - \M\\ ^ 1 the function f KM in Eq. ^ is 
reduced to the bare Tkm — Fkm, given by 



M 



C K CM 



x de 



fe - /, 



(55) 



S+U) 



see Eq. ([CI]). At high T > T and K £ M, the distribu- 
tion function can be considered smooth on the scale T; 



Eq. ( 55 ) reduces to 



Tkm(u,h,T) 



4Tk 
3w 



(fx- f M )F KM (6uj/r), (56) 
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FIG. 5. The dynamic conductivity according to Eq. (51 1 for 
different chemical potentials and T = 0. The cyclotron reso- 
nances are indicated by the solid arrows and the participating 
LLs, where L n denotes the n-th LL. The dashed arrows mark 
the disorder-induced transitions. The spectrum is calculated 
for a disorder strength a = 0.01. 



where we denote fx = I(Ek) and Suj = uj — Em + Ek- 
The function Fk m is given by Eq. ( C3 1 , see also Fig. [6] 
For T > T and K = M (intra-LL transitions), 



Fkk{u,H,T) 



IV F KK (8u/T K ) 



(57) 



3Tcosh 2 [(E K - E N )/2T] ' 
where N is the LL closest to the chemical potential /1. In 



the case of disorder-induced transitions, the conductivity 
vanishes in the clean limit T — > 0. 




FIG. 6. Solid line: The function Fkm, Eq. (|C3j), describing 
the shape of the disorder-induced conductivity peaks (\\K\ — 
\M\\ / 1) in the regime T,u c > F, see Eqs. (56} and (51 1. 
Dashed line: The function Fkm, Eq. (C6l, replacing Fkm in 
the case of the cyclotron resonance ( 



K 



I Ml 



strength and the width of the cyclotron resonance is reduced 
due to the vertex corrections as compared to Fkm that would 
result from the calculation with the bare polarization bubble. 
Here we assume K, M 7^ 0, in which case Fkm and Fkm do 
not depend on K and M. 



For the cyclotron resonance ||JV| — |M|| = 1, the vertex 
corrections are important leading to a more complicated 
e-dependence of the kernel; the corresponding functions 
T and F entering Eq. (51) are given in App. [cj The 
behavior of the functions F and F is illustrated in Fig. |6] 
The figure shows that the vertex corrections reduce the 
strength and width of the cyclotron resonances. 

We now discuss several limiting cases which describe 
the rich pattern of resonances in Fig. [5j 



A. Intra-Landau level transitions 

We start with the case u < T enabling only intra-LL 
transitions. This includes the dc limit uj —> and two 
distinct temperature regimes: a) T « T and b) T T. 
The contribution of the intra-LL transitions Lo —> Lo 
and Li — > L\ is illustrated in Fig. [5| 

a) When the temperature is the smallest scale, T <C 
r, only the level Ln, determined by the position of the 
chemical potential, contributes. The conductivity (51) 
becomes (ctq = e 2 /47t 2 ) 



V = CO J - -^NN(^, A*, J- ) ■ 



(58) 



In the dc limit uj — > 0, Eq. (58) acquires the forrrPH 
,2 



a = (T 



C N UJt 



UJ 



N 



1 



(£ - E N y 
r 2 

1 N 



(59) 
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Note that T is close to unity for small frequencies, hence 
the dependence on the disorder strength drops out for 
/i = Eff. In this regime, the dependence a oc r^ 1 oc 
a characteristic for classically strong magnetic fields, 
uj t\tt t r 3> 1, is exactly compensated by the increased DOS 
inside the LL: the average of vie) 2 over Lm is propor- 
tional to uj 2 /T 2 oca -1 . 

b) T > T, hence T > uj: With the help of the high T 
expression (57), the conductivity (51) reads 



^0 



E 



F nn (uj/T n ) 



- 3ulTcosh 2 [(E n -E N )/2T] 



(60) 



The summation limits are sent to infinity since the con- 
tribution of large energies where LLs overlap is exponen- 
tially small. 

In Eq. (60) the zeroth LL is special since its width 



is bigger by a factor of \J2 and its oscillator strength is 
enhanced by a factor of two. However, its contribution 
is significant only for N = and T < uj c . Note that 
for all other levels (n ^ 0) the function F nn (uj/T n ) does 
not depend on the LL index n. We obtain three high T 
regimes: 

b.l) For r « T « lun only Ln contributes, while the 
contribution from the levels farther away from the chem- 
ical potential is exponentially suppressed, 



N 



Suj^T 



(61) 



The conductivity a oc T 1 is proportional to the slope of 
the Fermi function in L N , f ~ 1/2- (e- E N )/4T. Note 



that the width of the peak in Eq. (61 1 for N = is bigger 
by a factor \[2 in view of r = Vzr. 

b.2) uj n < T < [i < e ov : As \i > T, the influence of 
zeroth LL can be neglected. Furthermore, at T > uj^ 



the sum in Eq. (60) can be converted into an integral, 
which gives 

64rl//l 3 

a = *o^^F NN (u>/T). (62) 

The conductivity in this regime is T-independent. 
b.3) For ujjy,fi<^T<e ov we obtain 

a = a 96 ^ 3 ) rr3 F NN (u)/T) , (63) 



where C{ z ) is the Riemann ^-function. We observe that 
the temperature takes the role of the chemical potential 



in Eq. (62). 



The dc-limit uj -S- of Eqs. (|61j), (|62j) and (|63J) is ob- 
tained using Fnn{0) — 1, and hence shows the same 
dependence on /i and T. 

The overall T-dependence of the dc conductivity 
(uj ->• 0) for T > T is shown in Fig. [7] (for [i ~ w c ) 
and in Fig. [8] (for fi ~ T). In both cases it shows a 
nonmonotonous temperature dependence. In Fig. [8] the 




3 

T/ui c 



FIG. 7. The temperature dependence of the conductivity for 
T 2> Tjv > uj and /i = uj c . The inset shows the small temper- 



ature behavior according to Eqs. ( 61 1 and ( 62 l (dashed lines) 



describing the 1/T decrease followed by a saturation into a 
T-independent regime. The thick line is calculated according 
to Eq. (60 1. The shaded area indicates the regime T < T, 



where the conductivity saturates at a /^-dependent value, see 



Eq. 1581 




FIG. 8. The temperature dependence of the conductivity for 
uj < r and \fi\ < To <C T. The dashed lines are the con- 



tribution from the zeroth LL [Eq. (61)] and all other levels 
[Eq. ( 63 1] . The thick line is calculated according to Eq. ( 60 1 . 



The shaded area indicates the regime T < F, where the con- 
ductivity saturates at a /i-dependent value, see Eq. (1581). 



T-indcpcndent regime ( 62 ) is not present since the cor- 
responding conditions cannot be met. At small T the 
contribution from L dominates. It decreases due to 
thermal smearing within the zeroth LL. With increas- 
ing T the influence of Lq decreases and the other LLs 
take over, which leads to an enhancement of the conduc- 
tivity due to thermal activation of higher energy states. 
In Figs. [7] and [8| the shaded areas indicate the crossover 
to the regime ( 58 ) , where the conductivity saturates at a 
/i-dependent value. The unusual T 3 and /i 3 dependence 
originates from the interplay between the energy depen- 



dence of the transition rates ( 53 1 and ( 54 ) and the level 
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spacing ([5| . It is therefore special for Dirac fermions and 
hence graphene. 

In the dc limit u — > 0, Eqs. ( 63 ) and ( 62 ) reproduce the 
results of Ref. where the regime T ^> ujn was studied 
for the dc conductivity. In addition, we find the 1/T 
behavior, Eq. (61 1, in the low-T range of T <C T <C ojm- 



B. Inter-Landau level transitions 

Figure [5] demonstrates that the nonequidistant LL spec- 
trum of graphene leads to a rich spectrum of resonances. 
Their strength depends strongly on the chemical po- 
tential. In what follows we discuss separately a) the 
cyclotron resonance transitions and b) disorder-induced 
transitions. 
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FIG. 9. The dynamic conductivity, Eq. (511, for fj, = and 
three different temperatures. The spectrum is calculated for 
a disorder strength a = 0.01. 

a) We start with the cyclotron resonances, \\K\ — 
\M\\ = 1. Apart from the transition Lk — > Lm also the 
transition L-m — > L-k needs to be taken into account, 
since it has the same transition energy. 

a.l) For T < T, Eq. (foil) reduces to 



/ x c K c M a uj ~ ~ 



if r M 



a.2) For T > T, Eq. (51 ) yields 



(64) 

Eq. ( [51) yields 

, v c k c m 4uj^ / Soj s 

0"(WJ = cr — {fK-jM+J-M-f-K)F K M I j; 

(65) 

Here fx are introduced below Eq. (56), and the functions 



T and F are given by Eqs. (C5) and (C6|. For M = 0, 
K = —1 the occupation of Lq drops out, 

a{uj) = a ° art ~ h) p -i>°( SuJ W ■ ( 66 ) 



In Eqs. (651 and (66), fj, and T enter via Jk,m only. 
The corresponding expression in brackets takes values 
between zero and one. 



b) \\K\ - \M\\ ^ 1: According to Eqs. (51]) and (|53j, 
the partial contribution of the disorder-induced transi- 
tions Lk — > Lm in the case K, M ^ is 

2(|iy| + |M|)[l + (|jfMAf|) 2 ] x , , rri\ 

<*KM = 00 1 , -rn9 -pj J-RMK^, H, J- )■ 



\{\K\-\M\Y-\f 
For M = or K = 0, 
c K c M r 2 



(67) 



CjfAf = CTO 



(K + 1)[M - (X + l)] 2 + (M + - (M + l)] 2 



\{M - K) 2 - If 



(68) 



where the function Tkm is given by Eq. (C2|. 

To make the result more transparent, we rewrite 
Eq. (67) for a particular set L-m — >• Lm (M > 0) of 



mirror transitions and for T»r, 

, > 4^ 2 r 



a w 



3w 2 f w 



(/— m ~ /m) F MM (6u/r) . (69) 



The function Fa/m is illustrated in Fig. Since 
2£ M , it holds 



M 



uT 

7^2 



(70) 



We see that the strength of the disorder-induced transi- 
tion is enhanced with increasing frequency ui oc yM. 

Our results are illustrated in Figs. [5] and [9] The 
strongest response corresponds to the cyclotron reso- 
nances. Indeed, Eqs. (64) - (66) show that the cyclotron 



peaks are a factor (JTJT 2 stronger than the disorder- 



induced peaks ( 67 1 - ( 69 1 . Additionally, we observe that 
the interband cyclotron resonance is suppressed by a fac- 
tor uj c /lj < 1 in comparison to the intraband cyclotron 
resonances, for which u c /lo > 1. 

Among the disorder-induced transitions, the intra-LL 
and the mirror transitions L—m — > Lm are the strongest. 
The latter become more pronounced with increasing fre- 
quency. The presence of the L\ —¥ L3 and the L\ — > L2 
peaks in Figs. [5jb) and (c) indicates that the chemical 
potential lies in the first LL. Furthermore, the intensity 
ratio of the L\ — > L 2 and the L — > L\ peaks provides 
information on the occupation of L\. For instance, the 
comparison of the calculated spectra in Fig.[5]to the mea- 
surements reported in Ref. 15 shows that in this par- 
ticular experiment the chemical potential was lying in 
L\\ since the Lq — > L\ resonance was stronger than the 
L\ — > Li resonance, we conclude that L\ was less than 
half-filled. No measurements so far have reported the 
disorder-induced transitions. However, a closer inspec- 
tion of Fig. 1 (a) and (b) from Ref. [T7] reveals a peak at 
roughly 140meV which should be attributed the disorder- 
induced transition L—\ — > L\ according to its position 
with respect to the cyclotron resonances L — > L\ and 
Li -> L 2 . 
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In Fig. [9] we illustrate the T-dependence of the con- 
ductivity for [i = 0. The strength of the cyclotron res- 
onance decreases to accommodate the increase of reso- 
nances at lower frequencies coming from thermal occupa- 
tion of higher LLs with smaller local cyclotron frequency. 
In particular, at moderate T = 0.3 lo c the L\ — > L2 tran- 
sition becomes visible, while for higher T it is absorbed 
into the background formed by other transitions. 



VI. SINGLE RESOLVED LANDAU LEVEL 

We finally consider the specific case of relatively dirty 
graphene with only one resolved LL at zero energy Lo 
within a quasi-continuum C. Due to the noncquidistant 
LL spectrum of graphene, this situation can still be re- 
alized for a moderately weak disorder strength, for in- 
stance, for a 



0.2, see Fig. 10 In the quasi-continuum, 
LLs are strongly broadened by disorder such that the 
quantization is effectively absent. Additionally, we con- 
sider probing frequencies much larger than the cyclotron 
frequency, lo ^> cj c , and assume U) ^> T, hence the tem- 
perature is effectively zero. 




FIG. 10. The DOS near the Dirac point in the sublattice 
b (sol id l ine) a nd a (dashed line) obtained numerically from 
Eqs. (loT and pll for a = 0.2. 



where vk=o{s) is given in Eq. (50) and uo{lo) — Lo/2nvp. 
The kernel describing the transitions C — > C is given by 



Kcc(£i,£2) 



aire 2 \e 1 e 2 \(el +e% - |eie 2 |) 



{e\ -e\± lo 2 ) 2 + a 2 Tr 2 (e 2 + e \ - | £l e 2 |) 2 ' 

which leads, up to small disorder-induced corrections, to 
the well-known universal value of the high-frequency con- 
ductivity in graphene, 



0-CC (w) 



1G 



57T 2 

16" 



a 2 Q{uj-2\e F \) + 0{a i ). (75) 



Multiplying the leading term in Eq. (75) by the de- 



generacy factor of 4, we obtain the clean conductivity 
a = e 2 /4h at B = 0, where we restored Planck's con- 
stant for convenience. Without vertex corrections, the 
conductivity acc obtains a contribution of the order a, 
which is canceled by the vertex corrections. The resulting 
disorder-induced correction to the universal conductivity 
of the clean graphene is of the order a 2 . 
The kernel ( 73 1 leads to the correction 



<7CoM = 



4 lo 2 ■ 

The total conductivity for lo > 2\ep\ is then given by 



(76) 



cr(w)l 



e 

16 



1 +4a 



+ 0{a 2 ) 



(77) 



It is dominated by the universal background with a small 
quantum correction oc B stemming from the resolved ze- 
roth LL. Despite its smallness, the correction oc B should 
be detectable in magnetooptical experiments since the 
background is £?-independent. The differential signal 



, der(u;) 



dLO 2 . 



UJ^>UJ C 



8lo 2 



(78) 



should provide an experimental access to the width of 
zeroth LL. 



We then distinguish two types of processes: (i) transi- 
tions C — > C within the continuum leading to a contribu- 
tion to the conductivity acc and (ii) transitions C — > Lq 
and Lo — > C giving the contribution aco- hi the case of 
lo > \ef \ +To, which means that we indeed probe Lq, the 
total conductivity <j{lo) = oco + °ec is given by 

f r ° de 
CTco(w)=/ K C o(e,e + Lo)— , (71) 
J-r w 

r- r ° He 
acc{u)= K cc (e,e + sign(|e| - \e F \)u)— . (72) 

Jr -u> w 

Here the kernel for the transitions C ^ L reads 

K C o(£,e + Lo) = — , (73) 



VII. SUMMARY AND CONCLUSIONS 

We have studied the linear transport properties of a 
single-layer graphene in the presence of diagonal white 
noise disorder and a moderately strong perpendicular 
magnetic field. Specifically, we obtained analytic results 
for the ac conductivity both in the semiclassical limit of 
high Landau levels (LLs) and in the quantum limit of well 
separated LLs. In both cases, several transport regimes 
are identified for different relations between temperature, 
the external frequency, the LL separation, and the posi- 
tion of the chemical potential. Additionally, we studied 
a specific case of a single resolved LL within a quasi- 
continuum of states nearly not affected by the magnetic 
field. Here, we show that corrections to the universal 
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value of the background interband optical conductivity 
in graphene should provide an experimental access to the 
effects of Landau quantization even in the limit of very 
weak magnetic field. 

The nontrivial topology of graphene, leading to the 
nonequidistant LL spectrum and to the unusual zeroth 
LL with strong sublattice asymmetry, makes its magne- 
totransport properties very different from conventional 
2D electronic systems with a parabolic spectrum. 

In the spectrum of graphene, both the quasiclassical 
regime of strongly overlapping LLs and the quantum 
regime of well-separated LLs can be realized simultane- 
ously. Which states dominate the transport properties is 
determined by the position of the chemical potential, the 
disorder strength and the external frequency. Two ex- 
treme limits are (i) the classical limit of strongly doped 
(or dirty) graphene where the effects of Landau quanti- 
zation are negligible and the conductivity is given by the 
Drude expression and (ii) the quantum limit where the 
dynamic conductivity is given by the set of delta-peaks 
with positions governed by the optical selection rules of 
the clean graphene. 

In the disorder-dominated quasiclassical regime the 
Landau quantization leads to quantum corrections to 
the semiclassical Drude conductivity. Since in high LLs 
JV > 1 the DOS is almost periodic with the period 
/VN <C w ci the Shubnikov-de Haas (SdH) oscil- 
lations in the ac conductivity <j(uS) are nearly harmonic at 
uj > <^, oc . At larger frequencies cj ~ cj c , the nonequidis- 
tant spectrum of LLs leads to the beating of the SdH 
oscillations which exhibit nodes for lo/lo c = \J1k + 1/2 
(k e N). Apart from the SdH oscillations, we also stud- 
ied quantum corrections that survive at high tempera- 
tures (T 3> To) where the SdH oscillations are strongly 
damped. These high-T oscillations also show a slow beat- 
ing at lu ^> uj 1 ° c . The nodes in the high-T quantum os- 
cillations are found to occur at oj/ui c = ykj2 (k G N). 

In the quantum limit of well-separated LLs, the con- 
ductivity is dominated by the cyclotron resonance (CR) 
transitions between LLs with indices N and M obeying 
\N\ - \M\ = ±1. For T <C uj 1 ° c , only a single transition 
Ljy — > Lm determines the conductivity. The spectrum 
of resonance peaks then strongly depends on the posi- 
tion of the chemical potential as seen in Fig. [5] Here, 
Fig. fBTb) corresponds to the measurement reported in 
Ref. Tl5l In particular, we can deduce the position of 
the chemical potential from the intensity ratio of the CR 
peaks. As for systems with a parabolic spectrum, the 
width of the individual resonance peaks for transitions 
between well-separated LLs, is given by 2(Tm + IV), 
(T > T), while the width of the Drude peak in the semi- 
classical regime is given by l/rt r . However, in graphene 
the ratio Tr tr = u c / ^J~a-K\e\ becomes energy dependent. 
The linear energy dependence of l/r tr has been con- 
firmed in measurements of the cyclotron resonance line 
width in graphene stacks. 17 Apart from the cyclotron res- 
onance, disorder enables transitions not respecting the 
clean selection rules in graphene. Their strength in com- 



parison to the CR is a direct measure for the disorder 
strength. The most prominent disorder-induced transi- 
tions are those between mirror symmetric LLs for which 
we find a oc u)T. The magnetoconductivity measure- 
ments reported in Ref. 1171 contain a resonance peak that 
we attribute to such a disorder-induced mirror transition, 
L_i — » Li (see Fig. [5]). For high T, multiple transitions 
can be addressed simultaneously. We studied the depen- 
dence of the low-frequency conductivity on the chemical 
potential and temperature and found that a oc /i 3 in 
the regime \i S> T and T 3> w' oc , whereas a oc T 3 for 
T 3> n,u) c . In contrast, the conductivity scales with 1/T 
for T <C 0J l ° c . The cubic dependence is a consequence of 
the Dirac LL spectrum. Specifically, it is a consequence 
of the energy-dependent transition rate which scales with 
/ [oj 1 ° c ] 2 in combination with the level spacing ui l ° c . 

Before closing the paper, we briefly discuss some of 
prospective research directions related to this work: 

(i) It woul d be very interesting to extend the experimen- 
tal studieal2I12 to the regime of overlapping LLs, where 
the beating of quantum oscillations should be observed, 
and to study systematically the dependence on temper- 
ature and doping discussed in Sec. [IV] The disorder- 
induced optical transitions in the regime of separated LLs 
also deserve a detailed experimental study. 

(ii) Theoretically, the effects of electron-electron interac- 
tions on transport and optical properties of graphene in 
moderate magnetic field, in particular, the interaction- 
induced damping of the magnetooscillations, is an inter- 
esting issue that has not been explored yet. Unlike con- 
ventional 2DES, in graphene interactions are generally 
strong and can directly affect the transport. Another pe- 
culiar property is that graphene supports fast directional 
thermalizatiorP^l due to a for ward scattering resonance in 
the electron collision integraP^U giving rise to effe ctive 
photon energy conversion via carrier mult iplicat ion JHSMQI 
This has tremendous consequences for potential applica- 
tions of graphene in rapidly developing field of graphene 
plasmonics and optoelectronicsPE^l 

(iii) The results of the present work, which studies the 
linear transport response of graphene to ac and dc per- 
turbations, serve a good starting point for studies of the 
nonequilibrium magnetotransport in strong ac and dc 
fields. First experimental work in this direction has re- 
ported the effects of heating of carriers by a strong dc 
bias on the SdH oscillations in grapheneP3l Just as dif- 
ferent cyclotron frequencies lead to a modulation of the 
linear-response dynamic conductivity, they are expected 
to manifest themselves out of equilibrium as well. Es- 
pecially near the Dirac point in graphene, the possibil- 
ity to create nontrivial nonequilibrium steady states at 
moderate B and high T, including population inversiorPH 
leading to optical gainj^l is tempting. From previous 
studies on semiconductor 2DES we know that nonequi- 
librium phenomena in moderate B are very sensitive to 
the details of disorder and interactions^ which are not 
accessible from standard measurements.^ Thus experi- 
mental studies of nonequilibrium magnetotransport phe- 
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nomena in graphene combined with their adequate theo- 
retical description should provide, in particular, valuable 
information about the nature of disorder and the role of 
interactions in electronic transport and relaxation which 
are still under debate. 



Appendix A: SCBA equation 

With the help of the Poisson summation formula the 
SCBA equation (11) is rewritten as the following equa- 
tion for e — E e : 



2 In 



-A 2 o A • e i27r l[e>?]'|sign([e?]'')A<: 

(e - Ef>) 2 + S1 § n ^ J > ' 1 _ A . e i2^|[ef]'| s i g n([ef]")/^ 



e(e-Sf) 



0. 

(Al) 



The coefficients of this quadratic equation also depend 
on the self-energies via 



ef = [ £ -e^)][ £ -e£( £ )], 



(A2) 



where /? £ {R, A}. The real part of <d £ is denoted as 
Q' e and the imaginary part as 6". We solve Eq. ( Al ) by 



iteration in A. Furthermore, for high energies we make 
the replacement 

-A 2 A 

In -f^ ^21n^±isign( £ )7T. (A3) 

(e - L, E ) \e\ 

Here ± refers to the retarded or advanced self-energy. To 



leading order in a the result is given in Eq. ( 23 1 



First we set all terms proportional to A equal to zero 
in Eq. (Al) and solve 



. a7rsign(e) 



(e-E e ) 2 [Z(e)Ti 



where the renormalization constant Z from Eq. ([15| has 
been used. Equation ( A4 ) yields 



£ ( e -s s ) + a ^ = o, 

(A4) 



e-E 



2r„ 



(A5) 



This solution is reinserted for into the oscillating part 
of Eq. ( Al ) and one obtains 



£-E = e+ ^-(l + 2^ A fe e i2,r 



fce 2 /"c 



(A6) 



k=l 



The leading terms in Eq. ( A6 ) yield Eq. (23 1 . 



Appendix B: Conductivity with vertex corrections 



1. Vertex corrections 

In this chapter we describe the SCBA in graphene and 
provide a detailed calculation of the conductivity in dis- 
ordered LLs of graphene including vertex corrections. 



(a) 



n, k 
n 2 , k 2 



-x-- 



m, k' 
m 2 , k' 2 



(b) 



Jni,mi 




n .ii 2 -in ■ "> i 



7712 





(d) Gf Gf Gf Gf 

[O-OMO-Oi 



U e+u> 



FIG. 11. (a) The disorder correlator in the Landau level basis, 
Eq. (B9 1, as a two particle operator, (b) The vertex correction 



r' 1 ) from Eq. (B12I. For brevity we omitted the momenta 



here, (c) The equation for the dressed verte x F . (d) The 
contributions to the conductivity kernel Eq. ( 35 1 including 



the vertex corrections denoted by the black triangle according 
to (c). 



We use the basis of eigenstates of the clean Hamilto- 
nian ([!]), 



IH-i) 



V2 \ s g n ( 



-0n=O,fe = <2 



_ Jky 





sgn(n)|0) 



(Bl) 
(B2) 



Here |n, k) are the harmonic oscillator eigenfunctions 
with origin l 2 B k. The matrix elements of the current op- 
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erator j x in the basis (Bl I, (B2) are 

(jx*)nk,mk f 



2L 



5kk' \pm$\m\,\n\-l + c n5\n\,\m\- 



(B3) 

where the coefficients c„^o = sign(n) and c„ = o = v2. 
With the self-energy from Eq. ([9| we obtain for the 
Green's functions in the case n, m ^ 

G n ,k;m,k' = $k,k' ^$njn[G +\ n ,k{£) + &n,-m [G-] n ,k (^) | > 

(B4) 

where 



[G+]n,k{e) = 
[G-]n,fc(e) = 



£ - (S„ + E fe )/2 + g„ 

(e-E )(e-E 6 )-£?2 
(E -E 6 )/2 



(B5) 
(B6) 



(e-£ )( e -E 6 )-.E» ' 
For n — 0, the matrix element of the Green's function is 

1 



Gn=o(£) 



(B7) 



e-E 6 (e) ■ 

The disorder correlator in graphene is defined as 

W(f - f) = (V(r) ® V(f)>dta , (B8) 

where V is the matrix of the impurity potential in the 
sublattice space. For diagonal white noise disorder we 
obtain for its Fourier transform the expression ©. The 
matrix elements of the disorder correlator in LL basis 
are formally obtained from the operator expression for 
the SCBA self-energy 0, 



W nk 



n2&2 \mk' ,ma k' 2 



(2tt)5 



(B9) 



X „(V'«,/c|e +I ' Z ' , '|V'n 2: fc 2 )^ CT (V>m 2 ,fc 2 |e "HVw), 



This equation is valid for a generic disorder in graphene. 
Here \tp n k}v denotes the ^-th spinor component of the 
eigenstates ( |B1[ ) and (B2|. For the vertex diagram in 
Fig. |TT|b) we exploit the relations between the LL in- 
dices in this diagram to evaluate the product of the phase 
factors (nk\ exp(iq • r)\n'k') in the case of white noise dis- 
order, Eq. ([8]). Using the solutions of the LL states in 
the Landau gauge for \n\ > |n 2 |, we find 



Jq 

(27T) 
Sk,k 2 



8ttP 



d< ^« e ip e (N-M+|m 2 |-|7i3|) / due~"u |riH ™ 21 



2tt 



f //.,| - I )\\ 2 ( (l TO 2|)!V L |n|-|n 2 |^ £ |n|-|n 2 |^ 



(\n\-l)\ 
(H)! 



(|m|)! 



> 2 |- 



n 2 \)\\ 2 ( (|m 2 |-l)! y |„| , na | h | | • 
(|m| - 1)! / I" 2 ' v ; I™ 2 !" 1 v ' 



(BIO) 



where L"(«) are the generalized Laguerre polynomials. 
Using the orthogonality of L"(it), we obtain the matrix 
elements for the diagram in Fig. [TT|b), 



W, 



nk,n 2 k 2 \mk' : m 2 k' 2 



SkM 8\n\-\m\,\n 2 \-\m 2 \ 



(Bll) 



x "^c n c n2 (5| n2 |^j_i_| m2 | -(- c m c m2 J| m2 |_]_i | n2 | 

The vertex function [Fig. [TTJb)] , to first order in the dis- 
order correlator, can be expressed as (/3,7 6 {R, A}) 

L n,k;m,k' ~ / j nk,n 2 k 2 ;mk' ,m 2 k' 2 



xG 7 , 

m 2 h 



kl,k\ (B12) 
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2 fc 2 ;mifei ( £ + W ) jm 1 k' 1 -n 1 k 1 G ni fcl ; „ 2fe2 ( £ ) , 

which gives 



+ c m 5\ m \ +1 .\ n \ IP p (e + uj,e) 



(B13) 



The ri-correlators are already for the case of overlapping 
LLs are given by Eq. (37). Their general definition is 



n^(ei,£ 2 ) = ^- ^ c m c m2 G^ n2 (ei)G7„ 2 _ m (e 2 ) 



n 2 ,m 2 



which can be written explicitly as 

2 oo 



(B14) 



2 ^o [ £ i - Es(ei)][ei - E&(ei)] - £^ +1 
e 2 -£ 7 ( £2 ) 



[e 2 -E2(e 2 )][ £2 -E 7 ( £2 )]-4 



(B15) 



For E a = £5 we recover Eq. (37 1. A closed expression 



for Eq. (B15) can be obtained in terms of the Digamma 
function W(z), which has been used in a similar context 
in Ref. \M 



u; c 2 [£i-^(£i)][£2-E2(£ 2 ) 



x {^{-^{e 1 )/^)-^{-^{e 2 )/u } l)-ut/@^^)} . 

(B16) 



Here G' 3 (e) is defined in Eq. \A2\ . Summing up all ladder 
diagrams, we obtain the dressed vertex 



Cm [n^( £ + ^, £ )]- 1 



(B17) 
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2. Conductivity bubble 



The conductivity kernel K(e,e + to) from Eq. ([35]) is 



the sum of all conductivity bubbles K^ 1 [see Fig. 11/1)], 
where /?, 7 € {i?, A}. At -B = 0, the dominant con- 
tribution away from the Dirac point and for intraband 
processes comes from the RA and AR bubbles while the 
RR and AA bubble give corrections of the order of the 
conductance quantum that are usually neglectedP^ At 
all four contributions are important. 
Without vertex corrections each bubble in Fig. [rT|d) 
is given by the expression 



0= £ £ 

ra,m,nj,TOj k,k' ,ki,k[ 

G 7 



jmk'.nk ^ n k,n2k 2 



(Bis) 



It is straightforward to rewrite this in terms of the II- 
correlators. 

= e 2 {n^(e,e + w)+n 7 ' 3 (e + w,e)} , (B19) 
which in the presence of vertex corrections turns into 

1 

(B20) 



O 



[II0T( E)e + w )]-l -mr 



1 



^(e + ^e)]- 1 -an 



Taking into account all diagrams from Fig. [njd) we ob- 
tain for the conductivity kernel, 



K( Sl ,s 2 ) = e 2 Re 



[n^(e 1 ,e 2 )]- 1 -an [II««(ei, e 2 )]-i - an 



(B21) 



a. Overlapping Landau levels 



In the regime of overlapping LL, the density-density 
correlators (B15l are calculated using the Poisson for- 
mula. Using the abbreviation 6f = (e — £^) 2 , we obtain 



+00 



n^( £l)£2 )=- [ £l -E^( £l )][ £2 -i] 7 ( £2 )] l dx 



+ 



k— — oo 



cxp(i27r&x) 



(01, - M - 1)(6 7 2 - |^|) 



£2-£ 7 (e 2 ) 



(B22) 



2tt [e 1 -^(s 1 )][(e 2 -zZ'>(e 2 )) 2 -L0 2 ] 



for the correlators (B15). Next, we neglect the second 



term in Eq. (B22) since it is small for energies e ^> u) c . 
The sum is dominated by the terms with k = and 



k = ±1 due to the presence of the coherenc e fac tor A , 
Eq. (|20|), in each term in the sum of Eq. 1B22I). The 



corresponding three integrals yield 



ttRR(RA) 

ei,£2 



2(ei-££)( ea -E 



R(A)x 



sign 



Im( £ i 



7- V" 



sign 



Im(e a -E«W) ! 



(B23) 



From the self-energy (|23[) we deduce that one may re- part 



place sign 



Im( £ 



— > sign( £ ) and we obtain the 



result (40 1 from the main text. 



Let us turn to the conductivity including the vertex 
corrections. Again, using the explicit form of the self- 
energies, Eq. (23), and expanding the correlators to lead- 



ing order in the coherence factor A, we obtain the Drude 



ReIL RA (e 1 ,e 2 ) = 



an\eie 2 \(el + £ 2 ) 



e 2 ±u 2 c ) 2 



a 2 n 2 {e\ + el) 2 



(B24) 
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The corresponding kernel ( 39 1 without vertex corrections 
reads 



i^ (0) (ei, 



£ 2) 



a-Ke 2 \e 1 e 2 \{e\ + e\) 



(£ 2 -£ 2 ±L0 2 ) 2 +a 2 ir 2 {£ 2 +£ 2 ) 2 ' 

If we include the vertex corrections we obtain 
1 



(B25) 



Re 



pI^( ei ,e a )]-i-o- 

a-K\exe 2 \{e\ + e\ 



|£i£2| 



(e\ -e\± ui 2 ) 2 + a 2 ir 2 {el + £ 2 2 - \£i£ 2 \) 2 



:V2 



(w ± w' oc ) 2 + 1/t, 



2 ' 
tr 



(B26) 



2r„ 



where we introduced the transport time r tr 
The factor 2 in the transport time can be understood 
within the Boltzmann theory from the Dirac factors 
(1+cos <fik,k') coming from the scattering matrix elements 
in addition to the transport factor (1 — cos 4>k,k') in the 
definition of the transport and quantum scattering time, 



i 



COS <pk,k' 

1 



\(k\V(q)\k')\ 2 . (B27) 



Here <f>k,k' is the angle between the initial and final di- 
rection of momentum. In the case of diagonal disorder, 
Eq. JB27} gives 



1/T ta 
1/7, 



1 — COS ( 

1 



v W(q) (1 + cos 4> ktk ,). 

(B28) 



From Eq. ( B26 1 , we obtain the result from Eq. ( 41 1 for 



the conductivity, which in the dc limit turns into 
1 Dr, 1 (fl^o)Dr tr 



tr(0) 



2tt 



2tt 



(B29) 



The higher order corrections in A yield results for the 
SdH oscillations and quantum corrections presented in 
the main text. 



b. Separated Landau levels 



Since LLs are separated and the kernel (B21) is pro- 



portional to the product of the DOS at energies e and 
e + to, the conductivity (34) can be written as 



KM 



de fe - fe+u 



K(s,e + Lu) 



(B30) 



x 6(r M -\e + u- E M \)Q(T K - \e - E K \) 



Comparing Eq. (B30) with Eq. (51), we identify 



•F K 



K 1 M 



de 



fe-fe 



x 6(r M -\e + u- E M \)Q(T K - \e - E K \ 
P(K,M) \ [U^ie^ + uj)}- 1 - air 



(B31) 



1 



p^e.e + w)]- 1 



where the factors P(K,M) are given by Eqs. (52)- (54). 
Without vertex corrections taken into account, Eq. (B31 ) 
reduces to 



KM 



uj 2 c k c m 



de 



fe - f 



€-\-LJ 



x 6(r M - |e + u - E M \)Q(T K - \e - E K \) 

Re {n RA {e, e + uj)- IL RR (e, e + u>) + {e o e + w}} 



P(K,M) 



(B32) 



leading to Eq. ( 55 ) in the main text 



With vertex corrections included, for the cyclotron res- 
onance \\K\ — \M\\ = 1 we obtain 



Re 



1 



1 



u 2 c {l B ^ 



~ an [I1 RR +U1 
2,T 2 v K {e) VM {e 



air 



(B33) 



9r 2 -. 



2S 2 )(S 2 - 2fc ) 



where we introduced the detunings from the center of the 
LLs Si = e — Ek and 5 2 = s + lu — Em and used T /ui c <C 1. 
The density of state s vk is given by Eqs. (49) and (50). 
The term Eq. (B31 ) with s and e + ui interchanged docs 



not contribute. For the cyclotron resonance \K\ — \M\ = 
— 1 the situation is reversed. The vertex corrections lead 
to an additional energy dependence in the denominator 
of Eq. ( |B33 ) in comparison to the expression 



Re II 



RA 

€,S + CJ 



n 



RR 



(B34) 

that would result from the bare bubbles in Fig. [TT|d) . 

Vertex corrections are important only at the cyclotron 
resonances (\\K\ — \M\\ = 1). In the presence of dis- 
order the selection rules of clean graphene are relaxed 
and the disorder-induced transitions (\\K\ — \M\\ ^ 1) 
become possible. For disorder-induced transitions, one 
needs to calc ulate the correlator Re (jl R f +UJ — H RR +U) ) 
entering Eq. (B32); the term with e and e + uj inter- 



changed is obtained by exchanging K and M. 

In the case \\K\ - \M\\ ^ 1, K,M ^ 0, calculation 
using T/uj c <C 1 gives 



Re(n^ +a) -n 



T 2 {e 2 



M 



RR \ 

S,€+U) ) 



vk{s)vm{£ + w) 



(4 - 4/ + ^c) 2 



(B35) 
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In the case K = M of intra-LL transitions, Eq. (B35) 
yields 

Re (n^ +w - n*? +a) ) «*(eW(e + «) 



(B36) 



For if = or M = and ||Jf| - |M|| ^ 1, 



Re (n*^, - n*» j 
. r 2 (4 + ^ 2 ) 



^(e)^M(e + w) 



2(4,-4+^)2 



(B37) 



Appendix C: Shape of the conductivity peaks in the 
regime of separated LLs 



shape of the conductivity peaks in the regime of sepa- 
rated LLs. 

For K 7^ M (inter-LL transitions) and T <C T, we 
obtain for Eq. ([B32l, 



r 

Fkm{u,ii) = - Q(T K + F M - \Slu\) 

UJ 

x Q(Sfj, + min{T K ,Suj + F M }) 
x 6(min{LA- — 6<j,Tm} — $V>) 

l.m.in{T K ,T M -Su,SiJ,}/T I 

x / deJl 

Jmax{-r K ,-rM-8u),6n-ui}/r V Cr 



x Wl 



(e + <wr) 2 

cm 



(CI) 



for the function (55). Here 5u) = w — Em + Ek and 
6(i = [i- E K . 



In the case K = M and for T < F, Eq. (B32 ) yields 




2 10 12 

Su/T 

FIG. 12. The evolution of the function 7bi(w, /Lt, T), Eq. |CT| ), 
for fixed T<T and for /x varying from n = — ro to /x = To 
inside the zeroth LL. 
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FIG. 13. The evolution of the function Jbi (w, /x, T) , Eq. |Cl| ), 
for fixed T < T and for /i varying from fj, = Ei — To to 
/x = _Ei + To inside the first LL. 



Using t he exp ressi ons (B33 )-(B37) for the Il-correlators 
in Eqs. (B31) and (B32), here we derive the explicit ex- 



Tkk{<»,h)= — e(2T K -\6u\) 

UJ 

x Q(Sn + mm{T K ,T K + Suj}) 
x 0(min{FA' — Sui, Yk} — <$/•*) 

x de \J\ - e 2 

Jmax{-l,-l-Su/T K ,5)j,/T K } 
Y2 



(C2) 



x y/l- (e + 5uj/T K y 
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FIG. 14. The evolution of the functions IFoi(u,[i,T) (top) 



from Eq. (CI I and ^01(0;, fi, T) (bottom) from Eq. ( C5 1 for 
fixed T <C r, both drawn for a chemical potential varying 
from /j, = to fj, — T/2. 



The functions entering Eq. ( 56 1 high T ^> T are 

'<5uW 



KM 



e(T K + T M -\5uj\) 



3r 

4T~ 



K 



r min{r K ,r M -fo}/r / 

< / deWl 

'max{-r A -,-r JU -5^}/r V °k 



(C3) 



pressions for the functions J 7 , F and J 7 , F describing the 



x Wl - 



(e + 5lo/T)i 
cm 
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In the case K — M (intra-LL transition) Eq. ( C3 ) turns 
into 



F, 



KK 



i-u>/r K 



de y/l -e 2 \ 1 - (e 



(C4) 



which can be expressed in terms of elliptic integrals, see 
e.g. Ref.HS] 
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FIG. 15. The evolution of the function J"n(o;, /Lt, T), Eq. | Clj , 
for fixed T < T and for /z varying from the value n = oj c to 
/i = uj c — 3r/4 inside the first LL. 



T«T, yields for Eq. (B31 1 



x 6(V + ndn{r K , 5u + F M }) 
x 6(min{rA- — 5oj,Tm} — <fy*) 

f .roin{r J f ) r M -5a!,5 i u}/r 3 A /l - 



de 



,{-iV,-r M -5u,,^- w }/r 9 + 8[e + 25w/r)][<iw/r - e] 

(C5) 



where an additional energy dependence enters in the de- 
nominator of the integrand in accord with Eq. (B33l 
This energy dependence also enters the function 



F, 



KM 



^)=e(r* + r M -|<M)^ 



min{r K ,r M -M/r 3a/1 



cm 



1 



de 



t{ -r K ,-rv-M/r 9 + 8[e + 25u;/r][fe/r - e] ' 

(C6) 



from Eqj57| for T > T. 



Figure. [12| illustrates the function Eq. (CI), for 
T<r. The arrow indicates the change with increasing 
chemical potential from /i = —To to fj, = +Tq. The func- 
tion J-qi (-^oi) shows an overall increase from (i = —T to 
// = +To followed by a decrease if the chemical potential 
jumps into the first LL and increases further, see Fig. 13 
Figure 
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compares the function Fqi, Eq. (C5|, to Fqi- 
We observe that the vertex corrections included in T re- 
duce the strength and width of the cyclotron resonance 
similar to functions Fkm aud Fk m illustrated in Fig. |6] 



Figure |15| shows the function Fu, Eq. (|C2j), for T < T 
and different chemical potentials. 
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